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Abstract 



We search and find no evidence for CP violation in r decays into the Kitv T 
final state. We provide limits on the imaginary part of the coupling constant 
A describing a relative contribution of the CP violating processes with respect 
to the Standard Model to be -0.172 < 9(A) < 0.067 at 90% C.L.. 
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The origin and source of CP violation in fundamental fermion interactions is a topic of 
great interest. CP violation has now been observed in the quark sector Increasing 
evidence for the existence of neutrino masses and their mixing opens the possibility of CP 
violation in the neutrino sector ||. It would be odd if the mixing effects were limited 
to the quarks and neutrinos only and did not appear in the charged lepton sector. Such 
mixing could lead to CP violation. There are strict limits on the mixing among the charged 
leptons coming from the searches for lepton number violation |J. Nevertheless, various 
extensions of the Standard Model allow for the existence of CP violation not only due to the 
mixing but also due to the interference between the W mediated and scalar boson mediated 
decays of the r to the same final state. This paper describes a search for CP violation 
in r decays. The results of the search are interpreted within the context of a model where 
CP is violated due to interference of IV-exchange and an exchange of a charged scalar with 
complex couplings. Our previous search for CP non-conservation in correlated r pair decay 
into 7T7r°z/ T final states || benefited from the large r — > 7i7i u T branching fraction yielding 
small statistical errors; however possible CP violating effects are isospin-suppressed for the 
mi u T final state JTO| , |n{] . In this search we study single r decays into the Kttv t final state. 
Although this decay mode has a smaller branching fraction, it is only supressed by the weaker 
SU(3) f symmetry and therefore has greater discovery potential. A previous search using this 



decay was reported in Ref. [12 



The most general way to search for CP violation is to define a CP-odd observable and 
then to determine its average value. A value different from zero would indicate CP violation. 
Various CP-odd observables have different sensitivity to CP violation. However, there is 
only one "optimal" observable £ that has the smallest associated statistical error [|13yi4f . 
For a decay described by CP-even P even and CP-odd P m components of the amplitude, the 
optimal variable is defined as 

£ -Podd/ -Peven- (1 ) 

In order to construct £ we need to know the explicit forms of CP-even and -odd parts of the 
amplitude in terms of experimentally measured parameters of the decay. This is possible 
only within a specific model. Thus the choice of £ is model dependent. 

We search for CP violation in the decay r — > Ktxv t in the context of a model where the 
CP symmetry is broken by an interference between Standard Model VT-exchange and an 
exchange of a scalar boson such as a charged Higgs |7||| with a complex coupling constant 
A. We assume that CP symmetry is conserved at the r pair production vertex. For this 
model, the matrix element for the r~ decay into Kix~v T final state is [[K|] 

A(t~ - Kn-u T ) ~ m(i/) 7m (1 - l 5 )u(r)f v Q» + Au(u)(l + <y s )u(r)f s M, (2) 

where fy and fs are the vector and the scalar form factors chosen to be Breit-Wigner shapes 
for if* (892) and Xq(1430) resonances, M = 1 GeV/c 2 is a constant providing a normalization 
of the scalar term, and is 

2 2 

= [(ptt-pkY 1 - fa +PK )2 b*+pxn ( 3 ) 

Here, p n , px, tu^, and are the momenta and masses of the outgoing pion and kaon. The 
square of the matrix element is 
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\A 2 \ ~ |/vf(2(g • Q)(Q ■ k) - (q ■ k)Q 2 ) + |A| 2 |/ 5 | 2 M 2 (g • k) + 2M(AMf s ft)Mm T (Q ■ k) 

- 2S(A)S(/ s /y)Mm T (Q-fc) , (4) 

where q and k are the 4- vectors of the r lepton and of the neutrino, respectively, and m T is 
the r lepton mass. The first three terms are CP even and the last, underlined term both 
violates SU(3) flavor symmetry and is CP odd. 

To construct the optimal observable we need to express (q ■ Q), (Q ■ k), Q 2 , and (q ■ k) 
in terms of experimentally measured decay parameters. From the energy and momentum 
conservation law we get 

(? • Q) = (Q ■ k) = -m^^f - ^][( < + K-^) 2 ) _ ^yft cos «, (5 ) 

Q 2 = 2m 2 n + 2< - [m% + (m 2 - m 2 K ) 2 ]/m 2 H , (6) 

(q-k) = (m 2 T -m 2 H )/2, (7) 

where m# is an invariant mass of the (jrK) system and a is an angle between the direction 
of a pion and the direction of a r lepton in a pion-kaon rest frame. The angle a is not 
measurable due to the unknown direction of the r. However, the cosine of this angle is 
statistically equal to the product of cosines of two other measurable angles in the pion- 
kaon rest frame: < cos a >=< cos f3 cos if) >. The brackets denote an averaging over the 
unobserved neutrino direction. The definitions of the angles (3 and if) in terms of measurable 
quantities can be found in Ref. [|i~0| , |Tll1 . 

In this study we use Eqs. (||)-(0) to express the CP-odd and -even parts of the squared 
matrix element in Eq. (f|). We use the CP-odd and -even parts of the squared matrix element 
in Eq. ([!]) to derive the optimal observable £. 

The data used in this analysis were collected with the CLEO detector at the Cornell 
Electron Storage Ring (CESR) operating on or near T(4S) resonance. The data correspond 
to a total integrated luminosity of 13.3 fb^ 1 and contain 12.2 million t + t~ pairs. Versions 
of the CLEO detector employed here are described in Refs. [IS] and [IB|. We estimate 



backgrounds by analyzing large samples of Monte Carlo events following the same procedures 
that are applied to the actual CLEO data. The generation of r pair production and decay 
is modeled by the KORALB event generator HT7|] , suitably modified to include the charged 
scalar contribution to the r — > K-rv t decay. The detector response is simulated with a 
GEANT-based @ Monte Carlo simulation. 

Tau leptons are produced in pairs in e + e~ collisions. Since the CLEO detector is more 
efficient for detecting Kg — ► 7r + 7r~ decays than for unambiguously identified kaons, we 
choose to make use of the r — ► Kgii ± h' T decay, which has a 3-prong topology. We select the 
candidate events on the basis of the one-vs-three topology. The one-prong 'tag' is based on 
a t candidate decaying into an electron, muon or a single charged hadronic track, and no 
more than one additional it . If one prong is identified as a lepton we require the presence of 
no more than one photon candidate with energy smaller than 100 MeV. The other, 'signal' 
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t is required to decay into a Kg, a charged pion, and a neutrino. We select events with 
four charged tracks and zero net charge. At CESR beam energies, the decay products of r + 
and t~ are well separated in the detector. Each event is divided into two hemispheres by 
requiring one charged track to be isolated by at least 90° from the other three tracks. Each 
track must have a momentum smaller than 0.85 -Ebeam to minimize background from Bhabha 
scattering and muon pair production. The momenta of all charged tracks are corrected for 
the energy loss in the beam pipe and tracking system. To ensure the existence of a K s decay 
in the three-prong hemisphere, the separation of the tracks in z at the r — <ft intersection 
must be smaller than 12 mm, the radial decay length of the K s candidate must be larger 
than 15 mm and the radial impact parameter of the K s must be less than 1 mm. 

Background from photon conversions is suppressed by requiring the cosine of the angle 
between two tracks to be smaller than 0.99. The invariant mass of tracks forming K s 
candidates must be within 12.5 MeV/c 2 of the known K s mass. To suppress background 
due to accidental combinations of the tracks we require the minimum impact parameter of 
the K s daughter tracks to be greater than 500 /zm, i.e., two times larger than the typical 
position resolution in the detector. To suppress background from r — > K s Kv T decay with 
the charged kaon misidentified as a pion we require dE/dx information for the charged track 
accompanying the K s to be consistent with a pion. 

To estimate backgrounds coming from r decays other than signal we use a Monte Carlo 
sample containing 39.6 million t + t~ events in which all combinations of r + and r~ decay 
modes are present, except for our signal process. Non-r background processes include an- 
nihilation into multi-hadronic final states, namely e + e~ — > qq (q = u, d, s, c quarks) and 
e + e~ — > T(4S f ) — > BB, as well as production of hadronic final states due to two-photon 
interactions. Backgrounds from the multi-hadronic physics are estimated using Monte Carlo 
samples which are slightly larger than the CLEO data and contain 42.6 million qq and 17.3 
million BB events, respectively. The background due to two-photon processes is estimated 
from Monte Carlo simulation of 37 556 27 — > r + r _ events, using the formalism of Budnev 
et al. [fT5|| . To study the CP- violating effects, we use Monte Carlo samples generated with 
and without CP violation [17|. The Standard Model Monte Carlo sample contains 170 000 
signal events which is four times that of the data, while CP- violating Monte Carlo samples 
generated with different values of the complex coupling A consist of 200 000 events each. 

To suppress the background from the multi-hadronic events (e + e~ —>■ qq) we require the 
invariant mass of the signal hemisphere to be less than the m T . To suppress background 
from two-photon interactions we require the missing mass scaled with the center of mass 
energy to be less than 0.65 and the scaled transverse momentum to be greater than 0.02. 
We also require the cosine of the angle between the beam-pipe and the direction of the 
missing momentum to be less than 0.95. Here, missing mass is the invariant mass of the 
difference between the 4-vector of e + e~ system and that for the total sum of all detected 
particles. Missing momentum is defined as a negative vector sum of all the momentum 
vectors of detected particles. The efficiency of the above selection criteria is (11.3 ± 0.1)%. 
A total of 11 970 events have been selected from the available CLEO data sample. 

After applying the above selection criteria to the Monte Carlo simulation of e + e~ — ► BB 
and e + e~ — > e + e~77 processes, we estimate that the remaining background from these 
sources contributes less than 0.2% to the data sample. The background from e + e _ — > qq is 
estimated to be (1.9 ± 0.2)%. The dominant background is due to misidentified r decays 
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with the largest contribution coming from the r — > KK°i/ T decay where we misidentify the 
charged kaon as a pion. We estimate this contribution to be (15.2 ± 1.7)%. The sources of 
background next in importance are due to r decays with a K$ in a final state with either 
a lost 7T°, such as r -> irKn°is T [(9.5 ± 1.0)%] and r -> KK°7r°u T [(3.1 ± 0.4)%], or a lost 
if£ in r -> K°7rKu T [(8.1 ± 1.4)%]. All other r decays contribute less than 1.0% each. 
The total background from r decays is estimated to be (39.2 ± 2.5)%, and from all sources, 
(41.3 ± 2.5)%. As a cross check of our signal selection procedure we calculate a branching 
fraction for r — > (Ktt)i = i/2V t an d obtain a value consistent with the Particle Data Group 
(PDG) |2(J within our statistical error. 

CP can be violated as a result of an interference between a vector [dominated by the 
If* (892)] and a scalar [e.g., the ^(1430)] resonance in the final state. To look for evidence 
of higher mass resonances we plot in Fig. [I] the invariant mass of the {Kgir) system for the 
data, signal Monte Carlo and backgrounds. We see no evidence for the if,*(1430) resonance. 
It can be seen in Fig. [I] that the K* mass peak in the data is shifted by approximately 
4.7 ± 0.9 MeV/c 2 with respect to the Monte Carlo simulation (which is based on the PDG 
mass of the K* + ). This is under study, but it does not affect the results presented in this 
paper. 

Another check for the CP- violating scalar component in the r decay is to look at the 
average value of the optimal observable as a function of the (Kgit) invariant mass. We 
expect the CP-violating effects to be maximal in the invariant mass range laying between 
the resonances, i.e., between 0.9 and 1.4 GeV/c 2 . In Fig. ^|we plot < £ > separately for r _ 
and r + as a function of the {K^tt) invariant mass for the data and for the Monte Carlo with 
maximum CP violation. A difference between the < £ > distributions for r~ and r + would 
indicate CP violation. We observe no difference in the < £ > distributions for the data and, 
therefore, no CP violation. 

To calculate the limit on the CP violation parameter A, we obtain the £ distribution for 
the data, for the Standard Model Monte Carlo simulation, and for the background Monte 
Carlo predictions. In Fig. [| we plot the £ distribution for both the full data sample and 
for the restricted region of the (Kir) invariant mass 0.85 GeV/c 2 < M(Ktt) < 1.45 GeV/c 2 
where the sensitivity to CP violation is maximal. Here, we change the sign of £ distribution 
for the r + decays to add r~ and r + samples together. The corresponding average values of 
< £ > are listed in Table |. 

TABLE I. Average value of the optimal observable in data, Standard Model and background 
Monte Carlo samples. 



Sample 


< £ >, 10" 3 


Full sample 


0.85 GeV/c 2 < M(Ktt) < 1.45 GeV/c 2 


data 

signal Monte Carlo 
r background Monte Carlo 
qq background Monte Carlo 


-1.5 ±1.5 
0.4 ±1.0 
0.6 ±1.6 
-18.1 ± 14.7 


-1.7 ±1.7 
0.5 ±1.1 
0.7 ±2.3 
-23.1 ± 19.1 


data (background subtracted) 


-2.0 ±1.8 


-2.3 ±1.9 
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FIG. 1. The {K° s -k) invariant mass for data (squares), signal Monte Carlo prediction 
(solid line) and background (hatched histogram). 
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FIG. 2. Average value of the optimal observable as a function of the {K^n) invariant 
mass for (a) data and (b) Monte Carlo with maximum CP violation 3(A) = 1. 




FIG. 3. The distribution of £ for the data (squares) compared to the sum of the back- 
ground (hatched histogram) and a Standard Model Monte Carlo prediction for the (a) whole 
data sample and (b) for the events with the mass of the (Kit) system ranging between 0.85 
and 1.45 GeV/c 2 . 
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To relate the observed mean value of the optimal observable < £ > to the CP-violating 
imaginary part of the coupling constant A, the 3(A) dependence of < £ > must be known. 
The optimal observable is a pure CP-odd quantity, therefore, its average value can be ex- 
panded in odd powers of the CP-odd part of 3(A). We have analyzed a Standard Model 
Monte Carlo sample with twice as many events as the data, finding an average value of £ 
consistent with zero (Table |). Therefore, the selection criteria do not introduce artificial 
CP violating asymmetry. A study of the CP- violating Monte Carlo sample showed that at 
small values of 3* (A) only the few first terms in the expansion contribute to the < £ > 

<£>~d3(A)+c 3 3(A) 3 . (8) 

We estimate c\ and C3 from the Monte Carlo generated with different values of 3(A) by 
fitting the observed average values of £ as a function of 3(A) to a cubic polynomial. The 
obtained coefficients C\ and C3 of Eq. (Bp are given in Table [TT]. 

TABLE II. Proportionality coefficients c\ and C3 for both full sample and the region of (Kir) 
invariant mass ranging from 0.85 GeV/c 2 to 1.45 GeV/c 2 . 



Coefficient 


Full sample 


0.85 GeV/c 2 < M(Ktt) < 1.45 GeV/c 2 


Cl 

C3 


0.0368 ± 0.0018 
-0.0135 ±0.0019 


0.0410 ± 0.0020 
-0.0127 ± 0.0022 



We use these empirically determined coefficients to estimate the value of 3(A). In Ta- 
ble |II| we list the results for both the full sample and the events with restricted (Kir) 
invariant mass range. 

TABLE III. The values of 9(A) as well as 90% confidence limits for both full sample and the 
region of (Kn) invariant mass ranging from 0.85 GeV/c 2 to 1.45 GeV/c 2 . 



Results 


Full sample 


0.85 GeV/c 2 < M(Ktt) < 1.45 GeV/c 2 


9(A) 

90% confidence limits 


-0.054 ± 0.049 
(-0.134, 0.027) 


-0.046 ± 0.044 
(-0.119, 0.027) 



To estimate the upper limit on the CP violating parameter 3(A) we first must estimate 
systematic errors. There are several possible sources of systematic errors that can contribute 
to this analysis. We treat these errors to be multiplicative if the source can modify the value 
of ci and to be additive if the source can bias the central value of < £ >. 

To construct the optimal observable we parameterize the scalar hadronic current as a 
product of a Breit-Wigner shape of the poorly known scalar resonance Kq(1A30) and a 
normalization constant M — 1 GeV/c 2 [Eq. (0)]. This assumption is a source of a systematic 
error on the calibration coefficient c\. We perform the following studies on the Monte Carlo 
simulation to estimate this error. We vary the width of the ^(1430) resonance by ~ 5cr 
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and re-calculate the coefficient C\. This results in a change of the value of C\ of ±4.4%. 
Similarly, we change the mass of the ^(1430) resonance from 1.35 GeV/c 2 to 1.45 GeV/c 2 
and obtain a variation of the value of c\ of ±11%. Thus, the overall conservative estimate of 
the systematic error due to the uncertainty in the mass and the width of i^o(1430) is ±12%. 

The choice of M = 1 GeV/c 2 is simple but not unique. Any Lorentz invariant quantity 
which has units of mass can be used as a normalization parameter. Another choice for M is 
the invariant mass of the (Kir) system. We generate a Monte Carlo sample with the (Kir) 
mass as a normalization parameter and re-estimate the value of C\ to be (0.0326 ± 0.0003). 
The resulting value of C\ differs from its nominal value by 2%. The overall multiplicative 
error on c\ due to a choice of the scalar current parameterization is ±12%. 

In Eq. (§) we define a Standard Model W exchange to have only a pure transverse 
vector current. However, chiral perturbation theory demands a scalar component |TIJ. It 



is assumed to be small and is neglected in TAUOLA. To estimate the systematic error due 
to the modeling of the W current we modify TAUOLA to include a scalar part. We use 
a Breit-Wigner shape for the 7^(1430) to describe the scalar component. The new value 
differs by ~ 3% from the nominal value. We take this as an estimate of the multiplicative 
systematic error due to the parameterization of the vector current. 

We estimate the values of the coefficients c\ and C3 in Table |T| from the Monte Carlo 
simulation. Therefore, the quality of the simulation may affect the result. We study the 
momentum distributions of the pion and of the reconstructed K s candidate in the signal 
r decay, both in data and in the Monte Carlo. We estimate a systematic error on the 
result by re-calculating the coefficients C\ and c 3 for the deviations between data and Monte 
Carlo parameterized as a slope of the ratio of the momenta distributions for the real and 
generated data. The multiplicative systematic error due to imperfect simulation of the data 
is estimated to be 9.3%. 

We study the systematic effects due to a possible difference of the track reconstruction 
efficiency for tt + and 7r~ as a function of the pion momentum. To estimate the size of this 
effect, we study the momentum distribution for charged pions in — > K s ir v T decay. The 
ratio of these distributions for r + and r~ decays is consistent with 1, and the maximum 
deviation characterized by a slope is fitted to be 0.01 ± 0.04. The introduction of such a 
slope to the data sample changes the value of 9(A) by ±0.009. We take this as a measure 
of an additive systematic error. 

A possible source of a bias is an asymmetry of the £ distribution induced by the remaining 
background. If we denote the number of signal and background events by S and B, then the 
contribution to the optimal observable due to the background is 

A <£>=<£> B B/(S + B). (9) 

Here, < £ >b is the value of the optimal observable in the background. We can estimate 
< £ >b from the r generic and multi-hadronic Monte Carlo simulations. For the r back- 
ground B/S is 0.413 ± 0.025 and < £ > B is estimated to be (0.6 ± 1.6) x 10~ 3 . Therefore, 
A < £ > is equal to (0.2 ± 0.5) x 10~ 3 . Such a change will modify the value of 9(A) by 
±0.014. Similarly, the background contribution from the multi-hadronic processes is esti- 
mated to be ±0.009. Therefore, the overall background contribution can modify the central 
value of 3(A) by ±0.017. 
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The overall multiplicative error is estimated to be ±15%, and the overall additive error on 
3(A) is ±0.019. Within our experimental precision we observe no significant asymmetry of 
the optimal observable and, therefore, no CP violation in r — * K-nv T decay. For a restricted 
range of the (Kir) mass (between 0.85 and 1.45 GeV/c 2 ) we obtain a value of the imaginary 
part of the scalar component in the r decays as 

3(A) = (-0.046 ± 0.044 ± 0.019)(1 ± 0.15). (10) 

The first error is statistical and the second is additive systematic. The overall expression is 
multiplied by the multiplicative systematic error. The corresponding limits are 

- 0.172 < 3(A) < 0.067, at 90% C.L.. (11) 

This limit is an order of magnitude more restrictive than that obtained in the previous 
search for CP violation in r — > Knv T decays. These results constrain the value of 3(A) 
at a comparable level to those from our study of t~t + — ► (7r _ 7r i/ r )(7r~7r i/ r ) M. However, 
the current result is about a factor of 10 more restrictive on the CP-violating parameters of 
Multi-Higgs-Doublet Models 0. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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